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Abstract 

A statistical-mechanical treatment of collision leads to a formal connection 
with transition-state theory, suggesting that collision theory and transition- 
state theory might be joined ultimately as a collision induced transition state 
theory. 

Reaction rate theories are divided into two approaches — collision theory and 
transition state theory. Differences between the two theories and attempts to either 
show their equivalence or to unite them into one have been a subject of discussions 
in the literatures and textbooks [1] . In this paper, I should like to give a discussion of 
possible joining of the two theories. Actually, whether implicit or explicit, we have 
met their relationship in unimolecular reactions. It is not my purpose to stress such 
a relationship, but rather to present a general theory which unifies the two theories 
and which may provide a framework of describing reaction rates. The discussion 
goes back to the very beginnings of the reaction rate theories. 

Consider a collision process between two molecules of A and B. We can discuss 
in the coordinate system of the center of mass the collision that occurs between 
the two molecules. All of the energy which goes into exciting the activated complex 
must come from the energy of relative motion of the reactants. Energy in the center 
of mass motion cannot contribute. According to the kinetic theory of collision, the 
rate constant has to be weighted by the Maxwell-Boltzmann distribution function 
f{u) of relative speed u, with integration over speeds from zero to infinity, to give 
the overall average rate constant: 
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fee = / o-uf{u) du, (1) 



where a is the collision cross section. The rate constant in this expression is given 
by M. Trautz in 1916 and by W. C. M. Lewis in 1918. 

It is instructive to evaluate the rate constant in terms of energy states instead 
of direct integration. We now consider the basic method of statistical mechanics of 
evaluating partition function [2]. Statistical mechanics states: 

The partition function is a sum over all states Q, very many of which 
have the same energy. One can perform the sum by first summing over 
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all the Q{E) states in the energy range between E and E + SE, and then 
summing over all such possible energy ranges. Thus 

Q = J2 e-'^"/'"'^ = J2 ^{E)e-^/'"'^. (2) 

n E 

The summand here is just proportional to the probability that the sys- 
tem has an energy between E and E + 5E. Since U,{E) increases very 
rapidly while exp{—E/kBT) decreases very rapidly with increasing E, 
the summand n{E) exp^—E/ksT) exhibits a very sharp maximum at 
some value E* of the energy. The mean value of the energy must then 
be equal to E*, and the summand is only appreciable in some narrow 
range AE* surrounding E* . The partition fimction must be equal to 
the value n{E*)eyip{—E* /ksT) of the summand at its maximum mul- 
tiplied by a number of the order of AE*/SE, this being the number of 
energy intervals SE contained in the range AE*. Thus 
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-^j, so lnQ = lnfi(E*)-^ + ln(^- 



'AE*\ , ^ , E* , , fAE*' 

6E 

(3) 

But, if the system has / degrees of freedom, the last term on the right is 
at most of the order of In / and is thus utterly negligible compared to the 
other terms which are of the order of /. Hence, the result agrees with 
the general definition S = ks lnO(£^*) for the entropy of a macroscopic 

system of mean energy E* . 

We have seen the basic method of statistical mechanics of evaluating the parti- 
tion function. If we apply this to the integration of Eq. (1), we would expect an 
expression for the rate constant to be 

kc=au*{^^y{E*)e-^'/'^-^, (4) 

where u* is a relative velocity for reaching the activated state. This summation 
indicates that the integration over the translational energy has a very sharp maxi- 
mum at the activation energy E* . The width AE* of the maximum, given by the 
square root of the dispersion, is very small relative to E* for a macroscopic system. 
The Maxwell-Boltzmann distribution function we have used is the one normalized 
to unity on integration over all states. For the results of more realistic calculation 
the normalization should be expressed in its explicit form. The following expression 
is then obtained: 



V 5E )\QaQb, 

This equation may equally be written in terms of an entropy change in reaching 
the activated state. As {AE* / 5E)Q,{E*) represents a number of energy states in 
the activated state, the expression for the rate constant may also be written in the 
form 

The basic method of statistical mechanics shows how an entropy term can be intro- 
duced in the kinetic theory expression. Here AS* represents the change in entropy 
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due to the change in energy in reaching the activated state. In a system of chemical 
reaction the entropy of the system is a function of energy E, vohime V, and the 
number of molecules N: S = S{E, V, N). Hence, we can replace AS* in Eq. (6), 
using the thermodynamic relations, by its generalization 

AS*^AS* — + t^-—, (7) 

for a system of chemical reaction. The rate constant can then be written in a general 
form 

kc = a«*e-^<^*/'=-^, (8) 

where AG* is the Gibbs energy change in going from the initial to the activated 
state. 

We can now uncover the collision theory expression for the rate constant of 
great interest. It has been stated that collision theory of reaction rates is not 
consistent with the fact that at equilibrium the ratio of rates in the forward and 
reverse directions is the equilibrium constant. However, it becomes evident that 
the kinetic theory of collision does not lack the entropy term that should appear 
in the expression for the equilibrium constant. Collision theory of reaction rates 
is consistent with the fact. Furthermore, it provides us with a kinetic theoretical 
derivation of thermodynamic expression for the equilibrium constant, by setting au* 
equal to that in the reverse direction according to the principle of detailed balance. 
We were able to show the thermodynamic equilibrium. 

The transition-state theory was published almost simultaneously by H. Eyring 
and by M. G. Evans and M. Polanyi in 1935. The rate equation for a bimolecular 
reaction derived by this theory of reaction is 



-E'/UbT 



(9) 
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from a partition function for the activated complex 

nv 

The partition functions Qa and Qb relate to the two reactants, and Q-^ is a special 
type of partition function for the activated complex. It is just like a partition 
function for a normal molecule, except that one of its vibrational degrees of freedom 
is in the act of passing over to the translation along the reaction coordinate. 

Equation (5) is very suggestive in relating collision theory to transition-state the- 
ory. The kinetic theory expression leads us to an idea of connecting with transition- 
state theory formula. By identifying with {AE* / 5E)Vt{E*), we can put both 
theories into some perspective. Prom their formal connection the reaction can be 
viewed as a succession of two steps — collision and transition state. The overall rate 
is then given by the sum of two average lifetimes: rate = {k^^ + k^g)~^- The rate 
reads explicitly 



rate = 
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The transition-state theory is concerned with a motion in the activated state and 
gives no explanation for reaching the activated state. The transition-state theory 
itself cannot be a complete treatment of reaction rates. A chemical reaction starts 
with an external interaction for reaction to occur first. It is very natural to consider 
molecular collisions as for an external interaction. The simple kinetic theory counts 
every sufficiently energetic collision as an effective one. Equation (10) suggests 
correcting the collision frequency by involving the translation along the reaction 
coordinate also in the evaluation of the transition over translational energy states. 
The essential feature of the argument is that transition state is brought about by 
energetic collisions and that the rate of a reaction is determined by the frequency 
of these collisions and by the resulting translations along the reaction coordinate. 

Indeed, the present treatment of reaction rates reflects the most important as- 
pects of unimolecular reactions [3]. Equation (10) is in exact agreement in form 
with the rate equation given by Rice-Ramsperger-Kasscl-Marcus (RRKM) theory 
of unimolecular reactions. The distribution function that has been used in RRKM 
theory is equal in expression to that given by the basic method of statistical me- 
chanics of evaluating partition function. But the present treatment has shown the 
rate equation in a general formulation of bimolecular reactions, and thus has given it 
a much wider applicability. The formalism provides a framework in terms of which 
molecular reactions can be understood in a qualitative way. The kinetic theory 
values are too high for all except atom— atom reactions. Hence, the transition-state 
theory values can be regarded as exerting important control over the rates of molec- 
ular reactions. It might be due to the high-pressure limit that has led krs to much 
closer agreement with experiment. 
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Addendum: Green's Function Method of Diffusion 



Pick's diffusion equation suggests using Green's function associated with the 
given partial differential equation. Actually wc can find in a mathematical physics 
book [4] a discussion of the Green's function for diffusion. But the book has very 
general discussion of the Green's function method. We must get down to practical 
cases. In this Addendum, we shall work through simple examples to see how it 
goes. 

We consider hydrodynamic derivation on an intuitive basis of the diffusion equa- 
tion. Let p{v, t) be the mean number of molecules per unit volume located at time t 
near the position r. If we assume that the total number of molecules in a volume 
V is conserved, we have, on account of the constancy of number, 

dN ^ Idp IdV ^ 

— =0, or — f- + 7-— = 0. (A.l) 
dt ' pdt V dt ^ ' 

In hydrodynamics, the rate of dilatation of the fluid is equal to the divergence of 

velocity [5]. If V-v = 0, the size of the volume element remains constant with time, 
that is, the fluid is incompressible. Equation (A.l) is thus written 

l^ + V-v = 0. (A.2) 

p dt 

Using the convective time derivative, we get the continuity equation in the familiar 
form. As a result, we can see the correspondence: 

— = 0«^=4>-^+V-J = 0. (A.3) 
dt dt ^ ^ 

When we use the diffusion equation J = —D\/p, we get the partial differential 
equation, satisfled by p. In the simple case of a dilute fluid, it gives 

I - DV'p = 0. (A.4) 

If there is a source or sink of particles, particles are not conserved and (A.l) 
becomes inhomogeneous equation that has a source or sink term, 

dN 

— = -XN. (A.5) 

dt ^ ' 

Hence 

^-DV'p=-Xp. (A.6) 

The form of the diffusion equation suggests using the Green's function associated 
with the diffusion equation for its solution. We need to construct the corresponding 
time-dependent Green's function that satisfies 

-^-DV^G=-S{r-r')S{t-t'). (A.7) 
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The Green's function is 

e-(r-r')'/4D(t-t') 

To interpret (A. 7) let G be the density of molecules in a medium. Then the impulsive 
point source involves introducing a unit of molecules at r' at time t'. The Green's 
function then gives the density at a future time for any other point of the medium 
and thus describes the manner in which molecules diffuse away from its initial 
position. It is a function of only the relative distance and the relative time between 
source and observation point. As in the case of the scalar wave equation, particular 
integral of the inhomogeneous diffusion equation (A. 6) is 

p(r, t)= j G{r, t; r', t')Xp{r',t') d^r' dt' . (A.9) 

By differentiating with respect to (r, t), we obtain 

(^^ - DV^^ p{v, t) = (J-^- DV^^ J G{r, t; r', t/)Xp{r', t') d^r' dt' (A.IO) 

= j - DV^^G{v,t;v' ,t')\p{v' ,t')d\' dt' , 

where uses have been made of VG = — V'G and dG/dt = —dG/dt'. Formal 
interpretation of (A. 7) may be obtained by considering the change of differential 
variables. 

The diffusion equation differs in many qualitative aspects from the scalar wave 
equation, and of course the Green's function will exhibit these differences. The most 
important single feature is the asymmetry of the diffusion equation with respect to 
the time variable. Multiplying (A. 6) by G and (A. 7) by p, subtracting the two 
equations, and integrating over space and over time gives 

j (g^+p^^ d\'dt'+ j D{pV'^G-GV'^p) dVdt' = -J GXpdVdt'+p{r,t). 

(A.ll) 

We may apply the Green's theorem to the second of these integrals. In the case of 
the first, the time integration may be performed. Finally 

J Gpd^r' + J d(^G^ - p^^ dSdt' + J GXpd^r'dt' = p{r,t), (A. 12) 

where d/dn is the normal derivative at the surface S directed outward from inside 

the volume defined by d^r'. G is chosen so as to satisfy homogeneous boundary 
conditions corresponding to the boundary conditions satisfied by p. While the 
first term of (A. 12) includes the effects of the initial value of p, other two terms 
represent the familiar effects of volume sources and boundary conditions. The 
various inhomogeneous equations and the initial-value problem can be solved in 
terms of the Green's function which satisfies homogeneous boundary conditions 
and a causality condition. 
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The problem discussed thus far is that of the absorption of one substance by 
another through which it can diffuse. This can be regarded as a problem in diffusion 
in which some of the diffusing substance becomes immobilized as diffusion proceeds, 
or as a problem in chemical kinetics in which the rate of reaction depends on the 
rate of supply of one of the reactants by diffusion. Chemical reactions are often 
considerable dependent on the mobility of the reactants as well as on the kinetics of 
the reaction itself. They are practical examples of processes involving simultaneous 
diffusion and reaction of one sort or another [6]. 

An irreversible first-order reaction will be an illustration. The rate of change of 
the total number of reactants is 

^ = -kN. (A.13) 
dt ^ ' 

By writing in terms of concentration, we obtain 

1^ - Z)VV = -kp, (A.14) 

and we get the solution in the form of (A. 12), with the rate constant k replacing A. 
Near-exact form of (A.14) has been used in the problems of conduction of heat and 
of carrier concentration in semiconductor. 

For an irreversible second-order reaction, the rate equation is 

^ = -fcTV^. (A.15) 

dt ^ ' 

The procedure (A.l) — (A.4) gives for second-order kinetics 

^ - DV'^p = -kNp. (A.16) 

In form this is "intcgrodifferential" equation. But wc can reduce this to a differential 
equation. The total number A'^ is a function of time. It is not a function of time 
and coordinates. The time derivative of N in (A.15) is differentiation following the 
motion of reactants. Thus the solution of (A.15) may be used to give 

dp r.^2._ J No 



_ DV^p = -k ^ p. (A.17) 

dt ^ \l + NoktJ'^ ^ ' 

We may then expect the solution in the form of (A. 12), but now with kNo/{l+Nokt) 
replacing A. Basically what we have done is to use the solution of the reaction rate 
equation to construct a solution of inhomogencous diffusion equation. It shows 
how the general form of equations can be deduced by using the solutions of the 
corresponding kinetics in the resulting diffusion equation. 

Usually more than one reaction contributes to the rate of change of a given 
species [7] . Organic reactions occur almost by a chain of mechanism. But once the 
resulting rate equation is given, its formal solution is obtained in integral form by 
use of the Green's theorem. The Green's function method shows a mathematical 
power of solving the problem of differential equation in terms of the integral equation 
formulation. One can see at once that it includes explicitly terms that depend upon 
the boimdary conditions and transport properties of the rcactant molecules. Thus, 
it is likely to be much better than the reaction rate equation, which provides no 
way to understand the transport process in a system in which the chemical reaction 
occurs. 
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